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Abstract 

We study the possible penetration of a static magnetic field in an idealized sample of 
many layers supporting a two dimensional charged chiral quantum fluid, to see whether 
there is a kind of Meissner effect. This is a non standard problem since the quantum 
fluid is incompressible having a gap in its spectrum. We find that the system shows 
an intermediate behaviour between superconducting and non-superconducting fluids, the 
magnetic field being screened or not depending on its orientation relative to the layers. 



1 



1. Introduction 



The study of the physical properties of non-conventional quantum fluids appears to be in- 
teresting and promising, even if a precise and detailed comparison with a real existing sys- 
tem is not attempted. In particular, attention has been devoted to possible non-relativistic 
quantum fluids living in a two-dimensional space, as a framework of models which could be 
related to layered superconductors (see for instance references [1] , [2] , [3] in particular for 
the magnetic properties of realistic layered superconductors). A general approach consists 
in introducing the problem by means of an effective lagrangian, a la Landau-Ginzburg, 
representing an universality class which should summarize the relevant degrees of freedom 
of some underlying microscopic theory not explicitly specified. A particularly interesting 
class, which is peculiar of two space dimensions, involves somehow a Chern-Simons gauge 
field, breaking parity and time reversal invariance, describing a non-relativistic quantum 
fluid which makes a chiral distinction between some left- and right-handed behaviour (ref- 
erences [4], [5], [6], [7]). 

Here we consider a particular universality class of that kind characterized by the fact 
that the spectrum of the fluid (before considering it charged and coupled to an electromag- 
netic field) has a gap. This chiral quantum fluid was presented and extensively studied in 
ref [7] (it turned out that it is formally similar to an effective lagrangian used to describe 
the Fractional Quantum Hall effect [8], [9], [10], [11]. In our case however, the physical con- 
text, the interpretation and the range of the parameters are quite different). It is of course 
very interesting to study the behaviour of the fluid when it is charged and can carry electro- 
magnetic currents, in order to compare with the behaviour of superconductors. Therefore 
here we address ourselves more specifically to a rather crucial issue. Due to the fact that 
this quantum fluid has a gap in the spectrum and it is like an "incompressible" fluid, does 
it possess the property of screening a magnetic field like a superconductor (Meissner effect) 
does? Since there is a gap, the standard mechanism similar to the Higgs mechanism in 
field theory looks impossible. There is no massless Goldstone boson which can provide a 
"mass term" for the electromagnetic field, thus resulting in a finite penetration length. In 
fact, in our case, apart from the collective modes representing the overall uniform motion 
of the fluid, there are not the linearly dispersing compressional modes, which are found in 
mean field treatments of fluids related to fractional statistics, see in particular [12], [13], 
[14], [15], [10], (for a possible different treatment, see appendix A of [7] and [16]). 

One can make this question more precise for the peculiar quantum fluid we are in- 
vestigating. Indeed, the reason why there is a gap is due to the fact that it is coupled to 
a Chern-Simons gauge field A cs . This field is determined by the constraint that its field 
strength V A A cs is proportional to the density fluctuation bp of the fluid (see section 
2 for the formulation of the effective theory). Since the fluid is non-relativistic there is 
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conservation of the total number of particles, that is J Sp = 0, consistently with the fact 
that V A A cs is "exact" , i.e. it is a total derivative. It is because of the CS gauge field 
that the Goldstone theorem is evaded and that the U(l) phase invariance breaks down, 
but there is a gap in the spectrum. 

Next, when considering the coupling with an electromagnetic field, the sum (A cs + 
eA em ) of the two gauge fields will appear in the covariant derivative. As observed in 
particular in reference [17], it looks like there remains an unbroken U(l) symmetry, because 
of a simultaneous gauge transformation of both A cs and A em which cancels in the sum. 
A "mass" for the electromagnetic field, and thus the Meissner effect, is recovered in the 
effective theory of reference [17] by means of a term which separately breaks the gauge 
invariance for A cs , due to the condensation of the "spinon pair" component. In our 
effective theory, instead, we consider a one-component quantum fluid and therefore there 
is no such a term. The dynamics is therefore quite non standard. 

In practice, we can put the question in the following way. The term possibly giving 
a "mass" to the electromagnetic field, which in the standard case would be proportional 
to (eA em ) 2 p (where p is the fluid density), is now (A cs + eA em ) 2 p. Therefore, there is 
the possibility that no "mass" and no screening of the magnetic field is obtained, if the 
configuration 

A cs = -eA em (1.1) 

is energetically favourite. (To avoid possible confusion with the literature on the Hall effect, 
we remind that in references [8], [9], [10], [11], the relation (1.1) implements the "filling" 
relation, the CS field being proportional to the total density and the electromagnetic field 
representing the strong, fixed, space-independent magnetic field of the Hall effect. In our 
case instead, the CS field strength is proportional to the density fluctuation, and the space 
dependence of the magnetic field is not a priori given, but has to be determined by the 
variation of the hamiltonian) . 

We consider here the case of A em corresponding to a static magnetic field, and the 
problem of its penetration in the bulk of a sample made of a stack of many two-dimensional 
layers, where the quantum fluid lies (see figure la-b). In ref [7] we already discussed this 
problem, finding a different behaviour depending on the orientation of the magnetic field 
relative to the layers. Namely if the field is orthogonal to the layers' plane, case of figure 
la, the system behaves like a type II superconductor with a finite penetration length. In 
the case of figure 1 6, when the field is parallel to the layers, instead the penetration length 
grows with some fractional power of the sample size. Thus in this case there is no Meissner 
effect, strictly speaking, in agreement with general considerations. 

In the study done in reference [7] we have included in the energy computation also the 
electrostatic effects, arising from the fact that locally bp ^ 0. We imagine of course that 
the layer where the fluid lies provides a uniform background that neutralizes the charge, so 
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that the fluid is globally neutral. But assuming that the background charge cannot move, 
there are locally electrostatic effects where 5p^0. These effects arise in particular in the 
configuration of equation (1.1), due to the fact that bp is proportional to V A A cs . The 
resulting additional energy helps in disfavouring the configuration (1.1). 

However it can well occur that, in more realistic cases, the background-neutralizing- 
charge is not really fixed and that it can in in some way compensate the local excess of 
electrostatic charge so that the system remains locally neutral. Thus it is important to 
reconsider the problem. 

Here we discuss the central issue of the screening of the magnetic field, by assuming 
that there are no additional electrostatic terms in the energy. Therefore we examine 
whether there are intrinsic mechanisms which would energetically disfavour the penetration 
of the magnetic field. 

Our result is that, even in this case, the pattern remains qualitatively the same as 
said above, namely we find screening, i.e. finite penetration of the magnetic field, in the 
configuration of figure 1 a and infinite penetration in the configuration of figure 1 6, although 
with some important quantitative differences. 

We have summarized in section 3 the main physical reasons for these results, while 
section 2 summarizes the formulation of the effective theory. Sections 4 and 5 contain the 
more precise and quantitative analysis. 

In conclusion, we find that a chiral quantum fluid, of the universality class described in 
section 2, besides having other interesting properties discussed in [7], behaves with respect 
to a static magnetic field as a "quasi-superconductor" , that is in a way which is somewhat 
intermediate between superconducting and non-superconducting fluids. This could open 
new prospectives in the scenario of the layered quantum fluids and the possible realistic 
systems. 
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2. The effective theory. 



The theory which we discuss in this paper is described by the following non-relativistic 
effective lagrangian density in two space and one time dimensions: 

C = i<j>*d <j> - t^I^I 2 - — (M 2 (2-1) 
2m 1 1 m 

here (j)(x, t) is a non relativistic complex field which plays the role of order parameter, 
related to the density by 

p=\<P\ 2 . (2.2) 

We will comment on the possible values of the dimensionless constant g in the next section 
(in reference [7] we used g = ^). 

We assume, as it is proper for a non-relativistic theory, a fixed total number of particles 
N, that is we keep fixed N = J d 2 xp. 

D is the covariant derivative D = V — % A cs and A cs is a Chern-Simons gauge potential 
related to the matter density fluctuation by the constraint equation: 
-> - 2n 

V A A cs = —5p Sp = p-p po =< P > ■ (2.3) 

k is a dimensionless number which we consider to be of the order of few units. As a 
consequence of the conservation of the number of particles we have consistently: 

d 2 xV A A cs = . (2.4) 



It is seen that the small deformations have the spectrum 



E(p) = v £2 + _L(£4 + 16^2) (2.5) 



where S = ^p is the gap. The whole spectrum, including the vortex excitations (see also 
reference [18]), and other relevant properties, in particular the chiral features, have been 
discussed in ref [7]. 

Now we are going to study the possible properties of the system of screening an external 
magnetic field. Since this is essentially a three dimensional phenomenon we suppose to 
build up a multilayered bulk of many two dimensional thin films separated by a spacing d 
(see figure 1). 

We are interested in stationary situations, and therefore we look for time independent 
configurations described by the hamiltonian 

H - 1 - ieA ' m - ilc >l 2 + \^ + ^ + i^ 2 } ■ ^ 

Here B and E are the magnetic and electric field, respectively. 
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3. Qualitative description and summary of the 
results. 



Before going into computational details, let us discuss the essential points of the magnetic 
field penetration problem, assuming that there is no electrostatic field. In a y-independent 
configuration, like in figure 1, A cs = Ay S {x) and A em = A y m (x) being the only non- 
vanishing components of the gauge fields, the phase of <fi can be taken to be a constant, 
conventionally zero, otherwise it would contribute an additional positive energy (of course 
it will play a role instead in the vortex configuration of section 4.2). Let us look first at the 
linearized form of the equations coming from the variation of the hamiltonian (dropping 
higher derivatives terms which are irrelevant for the discussion of zero modes): 

- d 2 T A em + ^(A cs + eA em ) = 
md ' 

_ 3t ( A V d lA cs + I*- (A cs + eA™) = (3 ' 1} 
md \2n J md 

One sees immediately that there is in principle a zero mode, corresponding to the configu- 
ration of equation (1.1). Thus if (1.1) could be competitive with the standard configuration 
it would spoil the screening of the external magnetic field. However to understand its rel- 
evance, one has to take properly into account the boundary conditions and to see how 
the allowed modes of A cs can actually implement (1.1). In the following subsections, we 
will analyze, in a rather qualitative way, the configuration (1.1) in the two geometries of 
figure la-b. A more quantitative study is developed in sections 4 and 5. Of course af- 
ter discussing whether or not the configuration of equation (1.1) is energetically favourite 
throughout the whole sample, we still have to find the optimum configuration and describe 
its space dependence. We will discuss it in detail in the following sections using variational 
methods, to be more general than the linearized approximation. 



3.1. Magnetic field orthogonal to the layers' plane. 

Let us start discussing the case of figure la, namely that with B orthogonal to the layers' 
plane. We suppose that the external magnetic field points in the direction of positive 
2-axis, and penetrates in the bulk in the x direction (see figure la). We further assume 
that the matter distribution is uniform in the y and in the z directions and we can choose 
the electromagnetic and Chern-Simons gauge field pointing in the y direction: 

^=(0,^,0) 

A cs = (0,A CS ,0) l ' J 
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and, in the gauge V • A em = V • A Ub = 0, depending only on x. With these choices the 
hamiltonian (2.6) becomes (L XiVjZ being the sample's sizes in the various directions): 

H = [ dx t 1 (\n -112 , l_ Aem , aCS\*.\ , 1 r?2 , 1 ^2 , <7 



LyL Z 

In this way we have 



{2S5(lW + M" + ^ OT | a ") + 5 SS + ^ + ^ 4 ")'} ' (3 ' 3) 



A em = xB 5p = 7^dxA cs . (3.4) 



From (3.4), using (1.1), we get 



S P = -^B (3.5) 



on every point inside each layer. 

Of course, the system is overall electrically neutral, therefore the support over which the 
quantum fluid lies in each layer acts as the neutralizing background. If we suppose that 
the neutralizing background charges are fixed, so that the fluctuation of charged matter 
Sp cannot be locally neutralized, then we have an electric field, inside each layer, given by 
the equation 

V • E = (3-6) 

here ^ is the three dimensional matter density (we remind that since dp is uniform in the 
z direction, the z component of E is always zero, thus E is effectively two dimensional). 
Recalling equation (2.3), V AA CS = ^-Sp, we see that E and A cs are dual two dimensional 
vectors that is 

E i = 7 ^e ij Af s i = x,y. (3.7) 
This means that the electric field contribution in (3.3) can be rewritten 

^ 2 - ^° S ? ■ (3-8) 

Therefore also when (1.1) holds we recover, through this electrostatic term, a "mass term" 
for the electromagnetic field ~ (A em ) 2 . This case has been extensively studied in [7] 
with the conclusion that the quantum fluid behaves, in this configuration, as a type II 
superconductor . 

Here we are interested in the case when the background charge structure is not so 
rigid, and we allow the system to neutralize, in some way, the fluctuation Sp. Thus, we 
analyze the behaviour of the system dropping the E 2 term in the hamiltonian (3.3). We 
can take into account the energy which is spent by the system for neutralizing the charge 
fluctuations, while retaining the same form for the effective hamiltonian, by considering 
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values for the constant g of equation (3.3) rather larger than the ones previously consid- 
ered in reference [7] (there, considerations based on anyon mean field theory suggested 

We begin by observing that from equation (3.5) we see that the cancellation (1.1) cannot 
hold everywhere on the two dimensional space. In fact if this were the case we would have 

f dx 6 P = A_f d ,VAl cs = -^- / dxB = -p-L x B (3.9) 
J d 2nd J 2nd J 2nd v ' 

which being different from zero violates the conservation of the number of particles. 
Thus, there must be somewhere an additional missing density 5p M . Actually, it will be 
concentrated on the edge of the sample, otherwise A cs would have a jump and (1.1) would 
no longer hold afterwards: 

8p M = —BL x 8{x-L x ) . (3.10) 

Notice that 

which exactly compensate (3.9): f (5p+5p M )=0. 

We see that Sp M is very large since it is proportional to L x which is macroscopic. There- 
fore a very large energy comes, for instance, from the term in the hamiltonian which is 
proportional to f(5p M ) . Thus, we foresee that the configuration of equation (1.1) will 
be severely energetically disfavourite, and that the quantum fluid will essentially behave, 
in the configuration of figure la, as a standard superconductor. This is confirmed in the 
detailed analysis of section 4. (It can also be that Sp M , so to speak, disappears because 
the fluid undergoes locally a kind of phase transition. But if the fluid is stable this too 
would cost energy, and the conclusion would be the same) . 



3.2. Magnetic field parallel to the layers' plane. 

Let us now turn to the second case in which the external magnetic field is parallel to the 
layers, that is points in the x direction and penetrates the bulk in the z direction (see 
figure 16). Here we suppose uniformity of the matter distribution along y. Again we can 
choose the gauge fields pointing in the y direction as in (3.2). With these assumptions the 
gauge electromagnetic field is 

A em = -zB (3.12) 



and, if (1.1) holds, 



dp = ^d x A cs = . (3.13) 
In 
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Equation (3.13) holds everywhere but at the border of the sample. In fact, we take formally 
the usual boundary condition that A cs is zero at infinity (that is outside the sample) and 
from equations (1.1) and (3.13) it is constant inside each layer, i.e. independent of (x,y): 

A cs = -eA em (z) 6{x) 6(L X - x) . (3.14) 

9{x) being the usual step function =0 for x<0, =1 for x>0. Therefore we have: 

Sp = -^A em (z) [5(x) - 5(L X - x)} . (3.15) 

Notice that now not only 5p = inside each layer, but also J dx^- = 0, as the conservation 
of the total number of particles requires. 

Notice also that the total amount of fluid accumulated at each border of every layer is 
J dx^j = ±^zB which remains finite for L x ^oo. 

We see thus that Sp is not macroscopically large and therefore we expect that its contri- 
bution to the energy will not be large. 

In reference [7] we have analyzed this configuration including the electrostatic energy which 
arises when dp is not neutralized by the background. In this case the electrostatic energy 
comes from the attraction of the two opposite charges accumulated at the boundaries. 
This gives a relatively weak effect, because the two boundaries are far apart, and the 
configuration of equation (1.1) remains energetically favourite. 

We will see that is even more so here, when we assume, like in section 3.1, that the 
background neutralizes also locally dp, and accordingly we forget the electrostatic effects. 
In section 5. we analyze this case in some detail and confirm that the configuration (1.1) 
is indeed favourite and find a penetration length l z , in the z direction, to be of the order 
of 

l z ~(L x 5) 1/2 (3.16) 

where L x is the sample size in the direction of B (we assume L y >L x ) and 5 is the thickness 
of the border region where the charge is accumulated (we remind, from reference [7], that 
including the electrostatic effect one finds l z ~ (L :E (i 2 ) 1 / 3 x (logarithmic corrections)) . Thus, 
in this case it is true that the system behaves differently from an ordinary superconductor, 
where l z is finite for L x ^oo. Notice however that still 4 s — >0 for L x ^oo. 
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4. Screening of the magnetic field orthogonal 
to the layers. 



4.1. Edge penetration. 

Here we take as a starting point the hamiltonian (3.3) dropping, as said, the electrostatic 
term: 

£-,=h {^('^ |2 + ieA ' m + - 4CS|2 ") + \ e2 + IM ■ (4 ' 1) 

We imagine that the region where the magnetic field is different from zero is, in absence 
of the sample, the interval — s<x<L x . The total flux of the magnetic field is fixed thus: 

$ f 

— = (L x + s) ■ Bm = dxB = fixed . (4.2) 

Ly J 

— s 

Then we redefine the zero of the energy subtracting the constant quantity — ^ J dx B 2 M , so 
that we can rewrite the second term in (3.3) as: 

\jdx (S 2 -Bl I ) = \jdx{B- B M f . (4.3) 

Therefore the new hamiltonian is; 

Now we suppose that the sample is placed with an edge at the origin of the x coordinate 
and that its length in the x-direction is L x . We imagine that the sample is much smaller 
than the region where the magnetic field is different from zero, that is s^>L x . 
Here for simplicity we treat the penetration of the magnetic field as if it were uniform, 
rather than exponentially decaying, and we call l x the penetration length. Since we have 
fixed the total value of the magnetic flux we have: 

^- = (s + l x )B = (s + L X )B M . (4.5) 

Ly 

This leads to: 

B = S -±^B M ~ B M , (4.6) 

S + l x 

and to: L L I 

B — B M = — —^-B M — B M = x x B M . (4.7) 
\s + l x J s + l x 
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That is: 



^Jdx(B- B M f ^ \b 2 m {L x - l x ) . (4.8) 



To check the meaningfulness of what we are doing, let us consider the case of the standard 
superconductor (A cs = 0, p = p ) and see what happens to the screening. 
The hamiltonian becomes: 

H ' L j ix {B-B M )\^'j ix( A^f . (4.9) 



LyL Z 2 



o 



Notice that since A em = xB, we have for the second term in (4.9): 







2md 

that is, using (4.8): 

Minimizing (4.11) with respect to the penetration length l x we get the standard value 



e 2 p 

Po/d being the three dimensional mean density. So our assumptions make sense. 
Substituting back in equation (4.11) we get the value of the energy of the standard con- 
figuration: 

\B 2 M (L x -h x \ . (4.13) 



L y L z 2 m V 3 



4.1.1. Detailed discussion of the configuration A cs = — eA em . 

Coming back to the study of our non-standard quantum fluid, as we said there is a possi- 
bility of ruining the screening of the magnetic field by means of the cancellation (1.1). If 
this happens we have: 
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which is far from being zero. 

Therefore in order to have the conservation of the number of the particles of the fluid, some- 
where there must be some missing fluid density, that we will indicate as 5p M , accounting 
for the mismatch. 

Due to the conservation of the number of particles, it follows that f dxdp = 0, thus it 
is not possible that (1.1) holds strictly, otherwise f dx B oc f dx dp 7^ 0. So we consider 
a configuration where (1.1) holds as much as it is possible that is everywhere but in a 
small region, say, of thickness 6. Thus we suppose that Sp M is concentrated in a small, 
microscopic, region 5 around x = xo, (as we said in section 3.1, xq^L x ), that is: 

ke I _ ( w ~ x o) 2 

<W = —B M L x —^e (4.15) 

Z7T Oy/n 

in such a way that 

dx5p M = —B M L X . (4.16) 



2tt 



Let us estimate the various contribution to the energy in this configuration. 
First: 

\ dx *J d ^ %r = «^' (4 - 17) 

where (/> = y/po + 8p and thus d x <j>=\ d f PM ~\H^. 
Second: 

±£* dx (-^B M+SPM ) 2 = ^(§)V, (j=Jf - l) . (4,8) 
We get the total energy 

" ±*b u ± + ±(£\*# u lJ4=±-x\ . (4.19) 



L y L z 4md 2tt 5 2 md\2TT J 5 



L 

We see that the energy gets a contribution proportional to Bj^-f-. Therefore compared 



' M S 

to the energy of the standard configuration, equation (4.13), we see that the configuration 
implementing the cancellation as in equation (1.1) has an energy which is larger by a 
macroscopic factor. 

Thus, we can disregard the possibility that the configuration (1.1) holds true throughout 
the whole sample. 
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4.1.2. A variational analysis. 



We study now in some detail the penetration of the magnetic field with a variational 
approach of the full hamiltonian (4.4). This will allow us to go beyond the linearized 
approximation, and to take into account possible important non-linear effects. Let us 
make the following ansatz: 

B = B e~ Xx A em = _^o_ e -\x 

A 

B M = Bq 

2n f x 2tv 
A cs = — J dx'8p{x') = -a—poxe-^ ( 42 o) 

dp = ^-d x A cs = -ap (l - px)e~^ 

Z7T 

p = Sp + po = po 1 - a(l - px)e~^ x . 

Note that J dxSp = 0. The role of the parameter a is to leave free the value of 8p at the 
edge. The goodness of this ansatz can be checked directly on the standard case (4.9) from 
which we get back the correct value (4.12). 

The numerical analysis indicates that p ^> A. Let us assume it for displaying a somewhat 
simplified expression, verifying a posteriori that p A is indeed realized. One gets to the 
following expression: 

md rr 27rei?oPn 1 ( ^ 2 Pn 1 . o , A 9 e 2 p B 2 1 3mdB?: 1 

- H = , 9 T a + ITT^ + l~9Po + Woc(a) a 2 + ' 1 



L y L z hp 2 A ' \2k 2 p 3 ' 4p» ru ' ■ 4 A 3 4 A 

(4.21) 

Here 

oo 

, ^ 1 /" , (2 -x) 2 e" 2x 

c(a) = - dx — ^ r 4.22 

V ; 8 7 l-«(l-x)e- a; V ; 

o 

is a slowly varying function of a, which for a small tends to a number of order of unit. We 
will consider it as a constant. Minimizing this expression with respect to a we find 

a = -\qB q , (4.23) 
A 

where Q = 27T 2 p 2 + ^2^ +4k 2^ c ■ Substituting back in (4.21) and minimizing now with 
respect to A we find: 



«-i-VSV 1+ ^ 0,+ S°- (424) 
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We see from equation (4.24) that the penetration length approaches the standard one for 
large g, for which Q — > 0. We recall again that a large value of g is expected because it 
effectively represents the fact that the system must spend energy to remain electrically 
neutral when Sp^O. 

Indeed this results have been confirmed through a numerical minimization with respect 
to A and a of the hamiltonian which is obtained from (4.4) using the ansatz (4.20), with 
/ u = y / po (corresponding to a coherence length of ~ 10A), for various values of Bq. Already 
for g — 1.5 we get l x ~1.2 ■ 10 3 A, and for g = 10 we get l x ~9.8 • 10 2 A (the standard value 
equation (4.12) is / x ~8.5 ■ 10 2 A). We have also numerically verified that taking p, smaller 
increases the energy of the configuration, confirming that /it>A as stated above. In fact, 
for p^oo we see that Q^~{^t thus a— >0, and l x — >\f^^- We see thus from equation 
(4.21) that the value p — >oo formally corresponds to the minimal energy. We have taken 
p at its physically reasonable maximum value, that is p ~ y^po- 

One can check that a is indeed small for that value, for g ~ 10 and _E?o<^10 3 gauss one 
gets a<10~ 2 (in all these numerical computations we have taken m to be the mass of the 
electron, ~ 250A _1 , and po = 4 • 10~ 3 A~ 2 ). 

We see that the penetration length is independent of the value of So and that a is pro- 
portional to —B (that is, A cs has the sign opposite to A em , as if the system would like 
the configuration (1.1) as far as it is possible). 



4.2. Vortices. 



In this section we study the penetration of the magnetic field from magnetic vortices. The 
starting point will be the following hamiltonian: 



\2md 



V - ieA e 



%A 



'CS 



2 + \ B2 + i>>' 



(4.25) 



We look for solutions of the form: 



0) = f(r)e 



ind 



(4.26) 



n is integer and represent the vorticity. 
Substituting (4.26) in (4.25) we get: 



H 



' d 2 r{-^— 
\2md 



1 



(d r f) 2 + -(n-erAr-rA c s Yf 



+ 1 2* 2 + id^ 2 }- (4 - 27) 
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Finiteness of (4.27) requires: 

/(0) = =s> 5/0(0) = -po 
lim 5p = =^> lim / = yfpo 



71 

lim A e B m = — 
i — >oo er 

lim A$ s = . 

1 — >oo 



(4.28) 



We solve this hamiltonian in a variational way with the following ansatz, which satisfies 
(4.28): 



2n 



er^ m = n[ 1 - e" A r ) =t B = -<9 r (rv^ m ) = — Ye 
p r 2 e 



2-\ z r 



2„2 



rA CS = _^ pnr 2 e -M 2 , 2 



5p=^-d r (rA^) = - P0 (l-^)e-^ 
ztt r 

f 2 =Po~ Po(l - p 2 r 2 )e-^ 2 . 



(4.29) 



Notice that J °° d 2 r dp = 0, therefore there is no missing 5p M . Later on, we will compare 
this configuration with a configuration where A cs = —eA em and we will need 5p M ^ like 
in the discussion of section 4.1.1. 

Notice also that in this case, differently from the case treated in the previous section the 
value of dp at the origin is fixed to zero by the requirement of finite energy, see equation 
(4.27), so there is no a parameter. 

With this ansatz we get, supposing p^>\, to be later verified: 

Minimizing (4.30) with respect to A 2 one finds easily: 

A 2 = P± ■ (4-31) 
2md v ' 

Therefore we recover the penetration length we had found in [7] studying the case including 
the electrostatic interaction: 

= \ = j¥± ■ («2) 
A y e z p 

Notice that this penetration length is independent of g. 
Minimizing (4.30) with respect to p we find: 



15 



which corresponds to a coherence length: 



1 / 2kn2 , A n*S 

r v = - = J . (4.34) 

p y 12imp + kp g 

Notice that, as is reasonably expected, the radius r v of the vortex decreases with increasing 
g. All these results have been confirmed through a numerical minimization with respect 
to A and p of the exact hamiltonian with various values of g. Notice that the results for 
the vortices agrees fairly well with those for the edge penetration. 

In particular for g = 2 we have r„~8A and / :E ~1167A, whereas for (7 = 10 we have r„~6A 
and the same value of l x (corresponding to the value of equation (4.32)). Notice that these 
results confirm that /it>A as stated above. 

Substituting (4.32) and (4.34) in (4.30) we find the energy of the configuration of one 
vortex: 

£^(^- 2 '^- 2 +¥)- («»> 

4.2.1. The A cs = —eA ern configuration with rotational symmetry. 

We now discuss, like in section 4.1.1, the possible cancellation (1.1) in the case of a config- 
uration which has a rotational symmetry, and the flux of B is given, like for the vortices 
case discussed above. We will compare it with the standard vortex configuration of the 
previous section. Following the discussion done there, it is unavoidable a missing 5p M . We 
suppose to have a fixed value of the magnetic flux. 



Therefore 



(B) = J d 2 rB = ^N . (4.36) 



J d \ 6 p=^ J d 2 rV AA CS = ~ Jd 2 rB = -^-nR 2 B = -hN , (4.37) 

where Ro is the radius of the sample, supposed to be a disk. 

So to have conservation of the number of particles we need a 5p M of the form: 

kN 1 (r-Hp) a 

where the thickness 5 is supposed to be microscopic. Note that 

d 2 r5p M = kN. (4.39) 
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We take Sp M concentrated at a macroscopic distance from the vortex core, say at the edge 
of the sample, that is at r = R . Otherwise, if it were concentrated at some different 
point, say at r = ro, then A cs would have a jump of ^ at r = ro, due to the fact that 
d r (rA cs ) = ^-rSp, and would no longer cancel eA ern for r>ro, contrary to our hypothesis 
that the cancellation holds throughout the macroscopic size of the sample. 
Let us estimate the most relevant contributions to the energy ^- coming from the presence 
of 5p M . We write: 

(4.40) 



H 



H 



+ 



H 



We find first 



Second: 



H 



2md 



d 2 r\Dd>\ 2 ~ 



8md 



+ 



d 2 r 



H 



{d r 5p M f 



8pt 



L Z J 9 md J 



Third: 



H 

l7 



md 



d 2 rB 2 



k_ 
2^ 



iV 5 



1 kN 
8md~6 2r 



7T R 2 



27T 2 N 2 



(4.41) 



(4.42) 



(4.43) 



We can distinguish two cases. 

a) iV = 1 (or few units). Then the most relevant energy contribution due to 5p M comes 
from {^-) v since the other two pieces are suppressed, at least, by the factor with R 
macroscopic. This energy (for 5 2 ~ is less or equal to the free energy of the standard 
configuration (4.35). Therefore for a small flux, that is for a very small magnetic field, the 
configuration where A cs cancels eA em is possibly favourite. 

b) Now we suppose B macroscopic, in other words B is fixed in the macroscopic limit 
Ro — ■> oo. Therefore from equation (4.37) iV ~ BRq — > oo. In this case the configuration 
where A cs cancels eA em gets the most relevant energy from (^-) 2 , namely from the piece 
proportional to j^. Therefore, for this configuration 

(4.44) 



H 

l7 



9 
4md 



ke\' 2 Rl , . 
— B — " + less important 
2tv J o 



We have to compare it with the energy of the standard vortex configuration, equation 
(4.35) multiplied by N, that is: 

< 4 - 45 > 

where c is of the order of por 2 : i.e. a finite number. Since ^ — > oo, clearly the standard 
configuration, or also a configuration of many standard vortices, is energetically favourite. 
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5. Weak screening of the magnetic field parallel 
to the layers. 



In this chapter we put our attention to the configuration in which the field B is parallel to 
the layers' plane (see figure 1 b) and study the screening effects in absence of the electrostatic 
interaction, as discussed in the previous chapter. 

With considerations very similar to those made at the beginning of section 4.1 we get, in 
this configuration, to the hamiltonian: 



— — = / dxdz 



1 



2md 



\d x cj)\ 2 +Cj\d z 4)\ 



+ 



2md 



em, aCS 



eA em + A 



(5.1) 

here Cj is a constant accounting for the Josephson coupling between the layers. 

We study the possible screening starting from the configuration in which we have the 

cancellation (1.1) inside the sample, that is: 



A cs (x, z) = -eA em (z) [9{x)9{L x - x)] . 
This yields (see section 3.2): 

6p(x, z) = ^d x A cs = -^ em [S(x) - 5(L X - x)] . 



(5.2) 



(5.3) 



Here, as in chapter 4, we suppose that the fluid density is confined in a microscopic region 
of thickness 5 so we approximate the 5-functions with 



Therefore we have: 



6p(x,z) 



S(x) 



ke 1 
2^7^ 



1 



(5.4) 



e s 



2 
J? 



A em (z) 



(5.5) 



Keeping this configuration, let us assume that the magnetic field penetrates in the z 
direction within a length l z and let us estimate it. We begin by estimating the various 
contributions to the energy in this configuration disregarding the terms proportional to 



e 2& ~0. First: 
1 



2md 



dxdz \ d x 4>\' 



1 1 (key 1 



4mdp y^2n V 2yr / ^ 3 Jo 



dz A 



\ 2 



(5.6) 



where we used the fact that d x 6—\ d ^ 5p "> ~ t-L^Ja). 
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Second: 



Third: 



Fourth: 



^ y dxdz (B - B M f = X -B 2 M L X {L Z - l z ) , (5.9) 



2 

(here to get (5.9) we have used arguments similar to those that led to (4.8)). 
Now putting it all together, and using the fact that A em (z) = —zBM, we arrive at: 

H 1 1 fkeVlf 1 A„ 2 ^ C J 1 fke\ 2 l n2 , 



^ y mdy/2^\2nj 5 \4p 5 2 * J M 3 8mdp ^ \2n J S 

+ ^B 2 M L x (L z -l z ) . (5.10) 
Minimizing equation (5.10) with respect to l z we find: 

'•=V¥- <6 ' n) 

This means that when the magnetic field is parallel we do not have screening in the usual 
sense because l z — > oo for L x — > oo. But we have a "quasi screening" in the sense that 
^->0. 

Substituting this result back in (5.10), we find the energy of the configuration: 



H 



^ \B 2 M L X L Z - l -B 2 M {L x f /2 {^)- 1/2 ■ (5.12) 

Let us compare this result to the energy of the standard configuration. That is, let us take 
A cs = and </> = - v /po; then the hamiltonian becomes: 

^ = ^(J^ +£ *-S>' <6 ' 13) 

which is minimal for the standard value (4.12) of the penetration length. The value of this 
minimal energy is: 

■{l z - 2 -I z ). (5.14) 



H - V T 
Y~ — 2 M l] 



Comparing this equation with equation (5.12) we see that (5.14) is always greater than 
(5.12) for L x macroscopic. So the configuration (1.1), in the geometry as in fig lb is always 
favourite. 
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